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Abstract 

We consider the linear growth of matter perturbations in various dark energy 
(DE) models. We show the existence of a constraint valid at z = between the 
background and dark energy parameters and the matter perturbations growth 
parameters. For ACDM 7o = ^|o lies in a very narrow interval —0.0195 < 7q < 
—0.0157 for 0.2 < i^mfl ^ 0.35. Models with a constant equation of state inside 
General Relativity (GR) are characterized by a quasi-constant j'q, for i^m,o = 
0.3 for example we have j'q ~ —0.02 while 70 can have a nonnegligible variation. 
A smoothly varying equation of state inside GR does not produce either |7ol > 
0.02. A measurement of 7(2) on small redshifts could help discriminate between 
various DE models even if their 70 is close, a possibility interesting for DE 
models outside GR for which a significant 7q can be obtained. 
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1 Introduction 



There is growing observational evidence for the late-time accelerated expansion of 
our universe [2]. This radical departure from conventional decelerated expansion is 
certainly a major challenge to cosmology. This non standard expansion could be due 
to an exotic non clustered component yet to be determined with a sufficiently nega- 
tive pressure called Dark Energy (DE). By analogy all models trying to explain the 
accelerated expansion are called DE models but many models go now well beyond 
this simple picture. While the usual Friedmann equations in the presence of a cosmo- 
logical constant term A seem to be in good agreement with the data, it is clear that 
other models with a variable equation of state are allowed as well [2] . While a cosmo- 
logical constant universe is appealing because of its simplicity it nonetheless poses the 
problem of the magnitude of the cosmological constant A. This is the basic incentive 
to look for other models where DE has a variable equation of state. An additional 
incentive comes from the possibility to have phantom dark energy at low redshifts as 
this excludes the quintessence models, models with a minimally coupled scalar field 
inside GR [3]. It might also be that one should change the theory of gravity, as for 
example in scalar-tensor models [H [5] , and a lot of research has focused recently on 
other modified gravity models and higher dimensional models. Sometimes, the back- 
ground expansion going back to high redshifts is enough to rule out some models 
but typically this is not the case: models of a very different kind will be able to have a 
viable background expansion where the low redfshift expansion is in accordance with 
SNIa data. 

Depending on the gravity theory one is considering, the growth of the pertur- 
bations, even at the linear level, will be affected. Indeed, while distance luminosity 
measurements probe the cosmic expansion, matter perturbations probe in a indepen- 
dent way (see e.g. [S]) the gravity theory responsible for their growth (and of course 
also for the cosmic expansion). The growth rate of matter perturbations could be 
probed with three dimensional weak lensing surveys (see e.g. [9j). Hence two DE 
models based on different gravitation theories can give the same late-time accelerated 
expansion and still differ in the matter perturbations they produce [ID]. This fact 
could provide an additional important way to discriminate between various models 
(see e.g. [11]) and it is therefore important to characterize as accurately as possible 
the growth of matter perturbations which is the aim of the present work. 

2 Linear growth of perturbations 

Let us consider the dynamics of the linear matter perturbations. These perturbations 
satisfy a modified equation of the type 



The gravitational constant Gefr depends on the specific model under consideration 
and the corresponding modification of gravity. For example, as was shown in 0, for 
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scalar-tensor DE models we have 



An equation similar to ([H) is also found for example in DGP models. The physics 
behind it is a modification of Poisson's equation (see e.g. [6j) according to (we drop 
the subscript m) 

^ = -47r G p (5 ^ = -47r Geff P 5 . (3) 

Of course, more drastic modifications are possible as well. In particular more elabo- 
rate DE models can be considered that could further increase the degeneracy between 
models inside and outside GR (see e.g. [12]). It is convenient to introduce the quantity 
/ = 4t2-^. Then the linear perturbations obey the equation 



f + r- + Ui-'^)! = l^a„. (4) 

ax z \ ax J 2Gn,o 

with X = In a. Equation (j4]) reduces to eq.(B7) given in [13] for = 1. The 
quantity 6 is easily recovered using / as follows 



S{a) = 6i exp 



f{x')dx' 



(5) 



We see that f = p when 6 oc a^, in particular / — > 1 in ACDM for large z and / = 1 
in an Einstein-de Sitter universe. 

An important issue is to characterize departures on small redshifts for different 
models. It is well known that for in a ACDM universe one can write 

f^^l, (6) 

with 7 = constant ~ 0.6, an approach pioneered some time ago [T^ and generalised in 
|15j . The characterization of the growth of matter perturbations using a parametriza- 
tion of the form ([6]) has attracted a lot of interest in the hope to discriminate between 
DE models based on different gravity theories. 

Of course it is possible to write in full generality 

/ = Qmi^y-^^ ■ (7) 

Let us consider the quantity 7' = For many models it turns out that 

7(2) ^ 7o + 7o ^ < z < 0.5 . (8) 

As we will see later, this could have interesting observational consequences. 

We now derive a constraint which is valid in general for any 7 (-2). It is easy to 
obtain the following equation 

- (1 + ^) ln^]^ + + 1(1 + 3(27 - 1) «;eff) = I UL'^ , (9) 
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where = wde ^de- From ([2]), it is easy to derive the following equation 



7o = [In ^mfl] 



1 1-1 



-^m,0 - 3(70 - -) ^eflf.O + - 7^ ^m,0 " ^ 

Z Z Ltjv 



(10) 



Equation ffTOl) is further simplified in models for which = 1 to very high accuracy. 
An example where this is the case is provided by scalar-tensor DE models for which 
< ^ - 1 < 1.25 X 10"^ We then obtain 



Gn,o 



-^Zo - 3(70 - I) w,s,o + I ^^U" - I 
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This does not mean that equation (ITT!) cannot differentiate between different gravi- 
tation theories satisfying = 1 but rather that if it does so it is through the value 
of 7o. This value is of course affected by the function Ges{z). We will assume below 
= 1 to very high accuracy. As we see from (fTTI) . we have 7q = 7o(7o, ^m,o, uide,o) 
which is clearly equivalent to a constraint of the form 

/(7o, 7o, f2m,o, wde,o) = . (12) 

In this connection one should note that fitting functions of 7(2;) proposed in the 
literature, even though they give a satisfactory fit for f{z) in models satisfying some 
assumptions, generically will not satisfy the constraint (I12p . In contrast the constraint 
( TT^ does not depend on any assumption about w{z). For fixed flmfl, wdea, there 
will be a value 7o,cr for which 7q = 0. However we will have generically 70 7^ 7o,cr and 
therefore 7o 7^ 0. 

Very generally, in any model for which the parameters ^^,0 and woe (and hence 
wde,o) are given, one can compute numerically the function 7(2) from the linear 
growth of the matter perturbations. Using ( fT2|) it is then possible to obtain 7q. We 
will do this in the next Section for various models inside GR. 

Before considering specific DE models, it is possible to derive some general con- 
sequences from the constraint (fT2|) . Generically 7q will not vanish, it needs not even 
be small. Let us consider 7g in function of 70 for Qm,o and wcEfi fixed. As we can 
see from Figure la, the constraint fll2p implies in excellent approximation a linear 
relation as follows 

7^, ~ c + 6 (70 - 0.5) 6-3. (13) 

The coefficients c, b depend on the background parameters b = b{wDE,Q, ^m,o) (re- 
membering that we take = 1). The coefficient b decreases while c increases when 
Qmfi decreases from 0.35 to 0.20 (see Figure lb). In contrast, c increases from —0.19 
for Qmfl = 0.3, to —0.17 for flm,o = 0.2 

For Qmfl = 0.3 we have c = —0.19. We stress that relation f|T3|) will hold indepen- 
dently of any particular model and is a consequence of the constraint f|T2l) . 

Depending on the specific model under consideration, for given background pa- 
rameters ^mfl and WDEfl: 7o will take the value 7o(7o) corresponding to the value 70 
"realized" by the model. Generically we will have 7q 7^ 0. 
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Figure 1: a) The left panel shows the constraint ( fTTi) for flm,o = 0.3 and various values 
ofwDE,o- We have from top to bottom: wdea = —^-4:, —1.3, —1.2, —1, —0.8. For 
given flmfi and wde,o, the couple 70, Jq is on the corresponding line for any model 
while 7q will depend on the value 70 realized in a particular model, b) On the right 
panel the constraint (ITT]) is shown in function of ^lm,o- From top to bottom we have 
wde,o = —1-2, — 1, — 0.8. We see that the coefficient b defined in (fT3l) increases for 
increasing Qrn,o and decreasing wde fl- 
it is also seen from Figure la that a small variation of 70, for fixed parameters 
u!DEfl, ^mfl, can induce a non negligible variation of 7q in accordance with eq. f|T3|) . 
In particular the relative change in 7q can be very large. We will show below that 
for wde = constant, the 79 values are restricted to a very narrow range with 7q 
—0.02. Even when one consider a smoothly varying equation of state, we still have 
-0.02 < 7^ < 0.005 (see below) for 0.20 < fi^,o < 0.35. In other words a smooth 
change in the equation of state of DE is not able to produce 7q < —0.02 for viable 
cosmological parameters. Therefore, a measurement of 7q outside this range could 
be a characteristic signature of a DE model where gravity is modified. Moreover, a 
precise determination of 7q could help to better discriminate between various modified 
gravity models. 

When VLmfl = 0.3 we have h = 3.13 for woEfl = —1, while b becomes smaller for 
WDEfl > —1 and larger for WoEfl < —1 (phantom DE today). Hence for wdea < 
we get a larger variation A7Q for a given variation A70. When Qmfl decreases, so 
does the coefficient b however this decrease is rather small for relevant cosmological 
values. It would be most interesting to investigate whether a precise determination 
of 7o is observationally accessible. In view of ([8]) this means that one should measure 
precisely 7(2;) on < z < 0.5. Another aspect concerns the extraction of one, or 
both, of the parameters flmfl or woEfl- If we assume erroneously that 7q = 0, a large 
error can result in the determination of fim,o or wde.o from the knowledge of 70. This 
is illustrated in Figure 2a. 
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Figure 2: a) On the left panel, the blue line shows the degeneracy in the ^mfi, "WoEfi 
plane for 70 = 0.555 assuming 7q = 0. The red, resp. green, dashed lines correspond 
to 7^ = -0.02 (top) and 7^ = 0.02 (bottom), resp. 7^ = -0.05 (top) and 7^ = 0.05 
(bottom). Ignoring the true non vanishing value of 7q increases significantly the 
uncertainty on the couples Qrn,o, wdea- b) On the right panel it is seen that models 
with very close 70 can be discriminated if 7 is measured for < 2; < 0.5 assuming 7 
is linear on small z, as often is the case. The lower the values of 70, the easier it is to 
discriminate these models through the difference in their slope 7q. For illustration, 
we have assumed here an error of 1%. 

3 Some specific models 

We now turn our attention to specific models inside General Relativity where DE has 
a known equation of state. 



Because of its simplicity and of the recent data that seem to imply that viable DE 
models should not be too far from ACDM (see however [IE]), this model plays a central 
role. We find for ACDM 0.554 < 70 < 0.558 (see Figure 3b) and -0.0195 < 7o < 
—0.0157 for 0.2 < fim.o ^ 0.35. Hence 70 varies very little in function of Qmfi while 
7q is negative with |7q| < 0.02. An observation outside these values, in particular a 
positive value for 7q, or a large negative 79, would signal a departure from ACDM. 

3.2 Constant equation of state 

We consider now a constant equation of state which includes of course the ACDM 
model. For the conservative ranges 0.2 < Qrn,o < 0.35 and —1.5 < wde,o < —0.5, 
we find 0.542 < 70 < 0.583 and —0.021 < 7q < —0.013. However, as can be seen 
from Figure 3a, for fixed parameter flm,o, the value of 7q is practically constant with 
7q ^ —0.02 for different constant wde despite a non negligible variation of 70. To 
summarize, for constant wde, 7o lies in the restricted range —0.024 < 7q < 0.01 while 
it is practrically constant if fim,o is fixed. However, as emphasized above (see Figure 



3.1 ACDM 
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2a), even in that case neglecting the true (nonzero) value of 7q can induce a significant 
error in the determination of Qm,o or wdea from 70. Finally, it is interesting to note 
that for given f2m,o all these models have essentially the same 79 while the parameter 
7o can vary by about 4% (see Figure 3a). 




Figure 3: a) The lines in colour on the left panel are the same as in Figure 1. The 
black line gives the true value of 70 realised in models with wde = wde^ = constant 
and Qjn.o = 0.3. It is seen that all models with wde = constant shown here have 
practically the same non vanishing 7q, 7q ^ —0.02. Note that 70 increases when wde 
increases, b) On the right, 70 is displayed in function of ^rn,o for the ACDM model. 



3.3 Variable equation of state 

Our analysis can be repeated for DE with a variable equation of state. To be specific, 
we take a smoothly varying equation of state of the type [TTl [18] 

wde{z) = {-l + a)+p {1-x) = Wo + wi , (14) 

1 + z 

where x = —. The corresponding evolution of the DE energy density can be computed 
analytically and yields 

Pde{z) = poEfl (1 + z)3(°+^)e-3^T*^ . (15) 

The results are displayed in Figure 4 for models with a negligible VLde ioi z ^ 1. For 
example, if we fix wq = —1.2, we can compute the values of 70 and 7q in function 
of (3 = wi and ^^,0- We find 0.55 < 70 < 0.56 and -0.022 < 7o < 0.005 for 
0.20 <nrr,o< 0.35 and < /3 < 1. 0. To summarize, a smoothly varying equation of 
state does not seem able to generate |7q| > 0.02. 

^Note that slightly lower values for 70 can be obtained for less interesting models with substantial 
phantomncss in the asymptotic past 
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Figure 4: The parameters 70 (left) and 7g (right) are shown in function of (3 = Wi and 
flmfl for a model with variable equation of state parameter wde = —l-2+(3 Hence 
all the points on the two surfaces have wcEfi = —1-2. The results for wde = —1-2 
are recovered for (3 = 0. We note for the left figure that 7q = is obtained for some 
particular combinations P, flm,o- 

4 Summary and conclusions 

Considering the linear growth of matter perturbations in various models, we give a 
constraint at 2; = 0, eq.l fTT]) valid for all models, including modified gravity DE models 
that satisfy = 1- This constraint implies that the quantity 79 is completely fixed 
by the remaining parameters 70, wde,o and flmfl- For the models considered here 
inside GR, |7o| < 0.02. Interestingly for models inside GR with constant wde, 7o 
is quasi-constant with 79 ~ —0.02 as the variation of wde,o is compensated by a 
simultaneous variation of 70 (for given fim,o)- 

We have generically 7o 7^ and we emphasize that a significant 7q could help 
discriminate between models, even if their 70 values are close. We have illustrated 
this schematically on Figure 2b. This potential resolution improves as flm,o goes 
up and/or wdea goes down and could be important when dealing with DE models 
outside General Relativity. We will give elsewhere specific models where this is the 
case [19]. Generally, this approach could be very fruitful whenever 7(2;) is close to 
linear on small redshifts < z < 0.5 so that the slope is essentially given by 79. So 
we feel it would be useful to try to measure ■y{z) on small redshifts, and not just 
7o. Finally it is important to realize that neglecting a small but nonvanishing 7q can 
induce a large error on the parameters flm,o, ujdea that one could extract from the 
growth of matter perturbations. 
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